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ABSTRACT 

Recent observations of white dwarfs in globular clusters indicate that these stars may 
get a velocity kick during their time as giants. This velocity kick could originate 
naturally if the mass loss while on the asymptotic giant branch is slightly asymmetric. 
If white dwarfs get a kick comparable to the orbital velocity of the binary, the initial 
Runge-Lenz vector (eccentricity vector) of the orbit is damped to be replaced by a 
component pointing toward the cross product of the initial angular momentum and 
the force. The final eccentricity may be of order unity and if the kick is sufficiently 
large, the system may be disrupted. These results may have important ramifications 
for the evolution of binary stars and planetary systems. 
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1 INTRODUCTION 



ISpruid (|l998l ) proposed that white dwarfs can acquire 
their observed rotation rates from mild kicks generated by 
asymmetric winds toward the end of their time on th e 
asymptotic gian t bran ch (AGB) l|Vassiliadis fc Woodlll993 ). 
iFellhauer et al] (|2003l ) invoked these mild kicks to explain 
a putative dearth of white dwarfs in open clusters (e.g. 
IWeidemannI 1 19771 : iKalirai et all l200ll ). The expected sig- 
nature of white dwarf kicks has been observed in M4 and 
NGC 6397 (jPavis et al.ll2006h . They found that young white 
dwarfs are less centrally concentrated than either their pro- 
genitors near the top of the main sequence or older white 
dwarfs whose velocity distribution has had a chance to relax. 

If white-dwarf kicks can have dramatic effects on the 
distribution of stars in a globular cluster, perhaps they could 
also affect other bound systems such as binary stars. The 
case of binary stars is interesting for a second reason as they 
represent a different physical regime than either white-dwarf 
kicks in clusters or neutron-star kicks in binaries. Whereas 
the latter two effects are impulsive, the white dwarf kick 
accumulates over many orbital periods, so its secular effect 
is more subtle to calculate. 



the decreasing mass of the system. Specifically, the angular 
momentum of the orbit, the Runge-Lenz vector, and the 
product of the binding energy and the period of the orbit 
are adiabatic invariants. 

Including the effects of an asymmetric wind requires 
further analysis. The total kick imparted by th e wind is on 
the o rder of a few to ten kilometers per second (jPavis et al.l 
over many tens of thousands of years, so the force 
exerted by the kick is typically much smaller than the grav- 
itational forces in the binary that result in accelerations of 
tens of kilometers per second per year — the kick is a per- 
turbation. Even in the limit that the kick can not be treated 
perturbatively but it is still adiabatic, one can use the fact 
that the Kepler problem with a constant external force is 
tractable in parabolic coordinates to find that the angular 
momentum and the Runge-Lenz vector along the direction 
of the kick and the product of the total mass of the system 
and the semimajor axis are adiabatic invariants. However, 
these three quantities are insufficient to determine the final 
state of the binary, most importantly, its eccentricity. 



2 CALCULATIONS 

The fate and evolution of binaries with asymmetric mass 
loss is complicated. Specifically the mass loss turns on and 
off over many binary orbital periods so its influence on the 
binary orbit may be adiabatic (in contrast with mass loss 
during the formation of a neutron star). If the mass loss is 
symmetric the orbit increases in size in inverse proportion to 



2.1 Centre of mass motion 

It is useful to separate the force that imparts the kick into 
a force on the centre of mass and a torque about the centre 
of mass. Using centre of mass coordinates, the Lagrangian 
of the binary including the kick on the giant star is 

1 . r 1 2 GMa M2 \ 

L = -MR- + -,r- + + F. . (r + -r) (1) 
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where M is the total mass, fi is the reduced mass, R is the 
position of the centre of meiss and r = ri — r2. The force 
exerted on a single giant star by its asymmetric wind is 



Fi = hv^ 



\M\ 



(2) 



where b is a dimensionless vector characterising the asym- 
metry; the magnitude of b can range from zero for a sym- 
metric wind to unity for a wind that blows only along a 
single direction lik a rocket. Therefore, the acceleration of 
the star is given by 

dvi 1 , I I 

^ = ^b«_,|M| 

and the total velocity kick for the single star is 

MwD,l 



(3) 



single 



= b«„ 



iln- 



Ml 



(4) 



so the force on the single star can be written in terms of the 
resulting kick 



Fi 



AVs; 



\n(MwD,i/Mi) • • 
From Eq. ([T}, the equation of motion for R is 
Fi = (Ml + A.f2)R 
and 

Av 

Single 

It ^ Ml -I- M2 \n{MwD,i/Mi) I I 
so the total kick to the centre of mass is 



At;, 



(1) 

binary 



At;singlc(Ml) 

\n{MwD,i/Mi) 



In 



M, 



WD,1 



+ M2 



Ml -I- M2 



(5) 



(6) 



(7) 



(8) 



where the product of the velocity and asymmetry of the wind 
is assumed to be constant (the mass-loss rate may vary) 
as the mass of the giant star decreases and ADaingie(Mi) is 
the typical velocity kick imparted to a main-sequence star 
of mass Ml. As the secondary becomes a white dwarf, the 
binary receives a second kick 



Ai;single(M2 



In 



MwD,i + MwD,: 
MwD,l + M2 



(9) 



A,;'"' - 

b-n-y In {MwD,2/M2 

If the mass of the primary and that of the secondary 
are about equal (due to dynamical biasing, for example, 
iMcDonald fc Clarkel[l993l ). the direct sum of the two kicks 
would equal the kick received by a single star. However, the 
two kicks will generally not be aligned with each other, so 
they must be added in quadrature yielding a smaller com- 
bined kick for the binary. The total change of the velocity of 
the centre of mass of the binary is typically about 70%-80% 
of kick imparted to an individual star (using the initial-final 
mass relation of llben fc Renzinill 19831 ). 



2.2 Constant Masses 

To treat the evolution of the orbit, the influence of the 
wind on the orbital dynamics must be determined, especially 
whether or not it is adiabatic. Fortunately, for the reasons 
mentioned earlier, the dynamics of the orbit with the asym- 
metric wind are remarkably similar to the Kepler problem. 
I Redmond! (|l964h found a generalisation of the Runge-Lenz 
vector in the presence of an external force (F = Fi). Specif- 
ically, the generalised Runge-Lenz vector is given by 



Lxp (rxF)xr 



where M is the total mass of the binary, /i is the reduced 
mass of the binary, L is the orbital angular momentum, r 
is the relative position of the two stars and p is their rel- 
ative momentum. The first two terms give the Keplerian 
Runge-Lenz vector. In this system only the component of C 
parallel to the applied force is conserved. The vector evolves 
according to 



-LxF. 



2GA/^i2 ^^^^ 

The generalised Runge-Lenz vector is no longer perpendic- 
ular to the angular momentum, 

LC^'^ (12) 
and the relative position of the two bodies is constrained by 



C r 



r-l 



(13) 



where I is the semilatus rectum. With some rearrangement 
this yields the equation for an ellipse with a focus at the 
origin in polar coordinates 

r = = ^ (14) 

l-ICIcose* l + ecose ^ ' 

where 8 is the angle between r and C. This equation shows 
that the vector C has a length equal to the eccentricity of the 
orbit and points from the pericenter toward the apocentre 
of the orbit. The external force naturally exerts a torque on 
the orbit, so the angular momentum is not conserved. 



L = r X F. 



(15) 



However, because the external force is taken to be a per- 
turbation, both C and L can be taken to be constant over 
an orbital period (P) . When averaged over an orbit the only 
component of L that remains must be perpendicular to both 
F and C because the secondary spends the same amount of 
time between pericentre and apocentre as between apoc- 
entre and pericentre and consequently the component of r 
perpendicular to C vanishes. The remaining component of 
the torque is 

^ C r „ „ I 



-C X F = 



-f- 

e2 Vl 



C-C e2 Vl+ecose* 

yielding the torque averaged over an orbit 

1 



- 1 C X F. 



C X F 



CxF 



I 
I 

Pe2 



1 + e cos 9 

1 de 



1 dt 



1 + e cos ( 



Using the definition of the angular momentum 
gives 



CxF 



1 



(1 + e cos ( 



(16) 



(17) 



(18) 



(19) 



To calculate this integral, the substitution z — e converts 
it to a contour integral, so 



cos 9 - 



1 f 1 



.d9 = 



dz 



(20) 
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which yields 
(L) = CxF 



CxF 



CxF 



1 fif 



dz 



(1 + f + ^) 



z dz 



I 



[z-z+) [z - Z-) 



where the contour is the unit circle and 

z± = - {±^l-e?-l 



(21) 

(22) 
23) 

(24) 



so z = z+ is a pole within the contour. Using the Cauchy 
integral formula yields after some algebraic simplifications 



(l) = CxF 



TT 2 



(25) 



LP e2(l- £2)5/2 g 
Using the value of the semilatus rectum from Eq. (|13p and 



Kepler's third law, P = ^/GMJo^ /{2n) where a = //(l-e^) 
is the semimajor axis, yields 



(L) = CxF 



l5/2 



2 + e^ 



a3/2 (1 „ e2)5/2 



= ^CxF. (26) 



To estimate the timescale of the evolution of the orbital 
parameters, Eqs. and (|26|) can be combined to yield 



C=5. 



9 



4 GMfi- 



(CxF)xF, L=- 



(LxF)xF (27) 



if the masses of the stars and the force are taken to be con- 
stant. If the prefactor is constant, these are the equations 
for a harmonic oscillator with a frequency of zero along the 
direction of F and a frequency of 



3F 

U! = - — 

2 



a 
GM 



perpendicular to F. The evolutionary timescale is 

forbital 



_ 1 _ 2 

^ ^ ^ ^ 3 F 



GM 
a 



^^winc 



(28) 



(29) 



When the prefactor can be taken to be constant, Eq. (I27|) can 
be solved by splitting the Runge-Lenz vector into a compo- 
nent parallel and perpendicular to the applied force, yielding 



C|| — Constant, 

/" :jLo'^F 

Cx = Cx,o cos ci^t -I- y 1 — Bg— smcot. 

LoF 

The angular momentum has a similar solution 
L\\ = Constant, 



= Lx.ocostjt + 



Lo CoxF 



F 



siriLut. 



(30) 
(31) 

(32) 
(33) 



In both cases the component of the vector along the direction 
of the force is conserved, and the perpendicular components 
vary harmonically. 



2.3 Variable Masses 

In general the masses in the binary and the force are not 
constant, so Eq. (Illf) must be modified to yield 



C = 



2GAf^2 
dt 



-LxF - 



GAf/i2 



(rxFj xr 

2GMn 



Lxp — ^(rxF)xr 



(34) 



If the applied force is written as F = {M2/ M)hv^i^g^\M\ 
where b is a dimensionless vector characterising the asym- 
metry in the wind and the mass ratio gives the force in the 
centre of mass frame from Eq. [T] Eq. (|34p can be expanded 
in powers of M. If Wwind and b are taken to be functions of 
the mass of the star, the first and third terms are first order 
in M, and the second and fourth terms are second order in 
M. The orbit average of the third term vanishes because 
the orbit is closed to lowest order in M. This leaves only the 
first term or Eq. (Ill|l at lowest order in M. In taking the 
orbit average to obtain Eq. (I26p . it was assumed that L and 
C were constant in time; therefore Eq. (|26p is only accurate 
to first order in M, so it is consistent at lowest order to use 
Eq. dnj and (f26)). 

In general the prefactor is not constant; however the 
combination Ma is an adiabatic invariant because the equa- 
tions of motion of the binary are separable in cylindrical co- 
ordinates as discussed in iJT] Furthermore, the applied force 
is proportional to M, so Eq. (|27|l can be written with the 
total mass of the system, M, as the dependent variable. This 
yields 



d^ ^ 9 



Ma 



dM^ 



4M2Af2 G 



nd (Cxb)xb 



(35) 



and similarly for the angular momentum. This assumes that 
the value of M is constant in time. One can argue that the 
mass-loss rate depends on the mass of the star, so changes 
in the mass-loss rate would only appear at higher order. 
This equation may be rewritten by changing variables 

to 



3; = In 1 -I- 



Ml) 



(36) 



where star 1 is the star that loses mass and the variable x 
increases in time. This yields. 



dx^^^V^hh) dx 



the equation for a damped harmonic oscillator with a damp- 
ing parameter that decreases in time. The original Eq. |35] 
yields the solution 



Cx(Mi 



Ml 



+ C: 



Ml 



where 



(38) 



(39) 



and the undamped frequency of the oscillator is 
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3 / Ma AUsingle 

2 y GM| ^Jorbital 



(40) 



and the parallel component is constant as before. This fre- 
quency is dimensionless because the eccentricity is now writ- 
ten as a dimensionless function of the mass ratios rather than 
time. 

The initial conditions give 



3.1 Weak winds 

Although it might not be obvious from Eq. (|38|l , in the limit 
where the wind asymmetry (b) vanishes, the eccentricity 
vector is constant. Specifically, in the limit of /3o ^ 1/2, 
Eq. (|38j through ((4T|) become 



Ci(Mi 



Cx,o-/?o 1 



Ml 



Ma 



A/i,o / A/i,o 



Lo X b 

Lob 



(48) 



Ci,: 



It 



Mo + Min 



±/3o' 



y/1 - 4/32M2 J ' 

Mo / 1-Cg Lpxb 



Mi,o V 1 - 4/3g L06 

If /3o < 1/2 then the solution is overdamped and 
provides a clear definition. On the other hand, if /3o 
then the solution is underdamped and oscillatory as 



(41) 

Eq. [381 

> 1/2 



Cx(Mi 



MMi 
MoA/1,0 

Cx.o cos 

(3o Mo 



1/2 



/31n 



/3 Mi.o 
^In 



Ml Mo 
Afi,oM 

Lo xb 







+ 

Mq + A/fi 
2/3M2 



Ml Mo 
Mi,oAf 



so in the limit where the asymmetry of the wind vanishes, 
the eccentricity vector remains constant. In general the ec- 
centricity will either increase or decrease slightly, but if 
one considers that the wind will propel the giant star in 
a random direction compared to the original, orbital angu- 
lar momentum, the final eccentricity is on average the sum 
in quadrature of the initial eccentricity and an increment 
due to the wind, therefore, slightly larger than the initial 
eccentricity. Furthermore, because on average the two sum 
in quadrature, the change in the eccentricity decreases as 
the initial eccentricity increases, so it would be difficult for 
such a weak wind to unbind a binary. 



3.2 Planetary systems 

^ In the case of a planetary companion, the value of /3o is 
large and approximately equal to /3 and M ~ Mi which al- 



(42) 



where /3 = /3o\/l - l/(2/3o)' 



lows further approximations. An important quantity in this 
analysis is the product of the undamped frequency and the 
mass of the secondary 



/3oM2 = O.O8M0 



1 km s 



M 



1 Mo 1 AU 



1/2 



(49) 



3 RESULTS 

The form of the evolution generally depends on the value 
of f3o specifically whether it exceeds one-half, and the initial 
and final masses of the objects in the binary. It is useful to 
obtain a more accurate value of /3o for various situations. 
Specifically, the rocket equation, Eq. |4l gives 



which is small compared to the mass lost if the velocity kick 
is small. An expansion in /3oAf2 for A'h ^ Mi, yields 



C_l(Mi 



Cx,o-/3o 1 



Ml 



M2 



A/1,0 / A/1,0 



■A 



AWsinglc — feUwind lu 



Ml,0 
Ml,WE 



(43) 



where At;siiigie is the magnitude of the velocity kick that the 
star would have received if it were solitary so 



/3o 



3 AUsin 



2 In 



0.1 



Mi,o 
Ai'single 



Ma 
GMi 



M 



1/2 



Mo 



(44) 



(45) 



1 kms-i V 2M0 1 AU; \IMq 



Au 



single 



M 



1 km s^^ V 1 Mq 1 AU 
where Mj is the mass of Jupiter and 
Ml , Mi,o 



In 



Ml 



In 



O.38M0 + 0.15A'/i,o 



(47) 



was taken to be 0.63, the value appropria te for one solar 
mass using the initial-final mass relation of llben fc Renzinil 
l| 19831 ). 



which is identical to the weak wind case, Eq. (|48p even 
though the undamped frequency is large. In both cases, the 
important factor f3oM2 is small yielding a similar expansion. 
Specifically if the initial eccentricity is low, as one would ex- 
pect for planetary systems, the final eccentricity of the sys- 
tem will increase linearly with the strength of the wind and 
the amount of mass loss, resulting in final eccentricities of a 
few percent. In general the initial eccentricity and the incre- 
ment due to the asymmetric wind will add in quadrature, 
resulting in a modest gain in eccentricity due to the wind. 



3.3 Strong winds 

In the case of a strong wind (a large value of Ansingio com- 
pared to the typical orbital velocities), the value of P0M2 is 
large, as is the value of /3o- Again /3o is approximately equal 
to f3, but the two terms that multiply the sine function in 
Eq. (|4 2l) are similar in m agnitude. The results of lDavis et al.l 
(120061) and iHevll (|2007l ) indicate that the kicks that white 
dwarfs receive as giants could be as large as 5-10 km/s, so 
/3oM2 ~ IM0, so the approximations in i]3.1l and l3.2l do not 
apply. 

As before the final Runge-Lenz vector is a linear com- 
bination of the initial Runge-Lenz vector and a vector of 
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length about unity in the direction of the cross-product of 
the initial angular momentum and the asymmetry. In gen- 
eral the results depicted in Fig. [1] are rather complicated. 
The final Runge-Lenz vector is given by 

C|i =C||.o,Cx+^c^Cx,o + ^Loxb^^^^ (51) 

Loo 

where the two coefBcients are depicted in Fig. [T] 

The two contributions are about ninety degrees out of 
phase, so the final eccentricity may have little memory of its 
initial value. Furthermore, even if the initial eccentricity was 
small, the second term is of order unity, so the final eccen- 
tricity may be large. For a particular white dwarf binary one 
can infer the initial mass of the white dwarf from its current 
mass using, for example, Eq. (|47p and backtrack to obtain 
the initial semimajor axis and estimate the contribution to 
the final eccentricity from kicks of different magnitudes us- 
ing Eq. (|42l>. 

3.4 Disruption 

In Fig. [T] the length of the final eccentricity vector may 
exceed unity; this implies that the binary has become un- 
bound. To estimate the velocity kick required to unbind bi- 
naries with various initial properties the maximum eccen- 
tricity of an initially circular orbit is determined. For small 
velocity kicks, this typically occurs at the end of the mass 
loss, but for larger kicks the peak eccentricity according to 
Eq. (|42|l may be reached while the mass loss is ongoing. 
Fig. [2] depicts the velocity kick that the primary would re- 
ceive if it were single that is large enough to unbind the 
system. For all of the primary masses, the required velocity 
kick increases with decreasing mass of the secondary. This is 
simply due to the translation into the centre of mass frame 
of the system. For large primary masses, the primary loses 
approximately 85% of its mass independent of its mass, so 
the results become nearly scale -free in this regime, the re- 
quired velocity kick is about GMo/(2ao) for equal mass 
binaries and increases in inverse proportion to the mass of 
the secondary. The dependence becomes shallower when the 
final mass of the primary exceeds that of the secondary. 
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4 CONCLUSIONS 

The possibility of asymmetric winds in asymptotic giant 
stars opens a range of new phenomena in the evolution of 
binary stars. The presence of an asymmetry in the wind typi- 
cally increases the eccentricity of the binary, so a comparison 
of the properties of a binaries containing a white dwarf to 
binaries of main-sequence stars could provide an indepen- 
dent probe of the importance of white dwarf kicks. In equal 
mass binaries if the wind is sufficiently asymmetric to induce 
a kick on the order of the orbital velocity of a single star, 
the kick is likely to unbind the system, in contrast with the 
symmetric mass loss that only results in an increase in the 
binary separation, so the binary fraction of main-sequence 
stars and white dwarfs may differ providing an additional 
probe of this process. The effects of white dwarf kicks in 
planetary systems are more subtle. Even large kicks do not 
typically unbind the system; however, kicks may pump the 
eccentricity of the orbits to high enough values to change 



_ I I I \ \ \ I \ \ \ \ I \ \ \ I \ \ \ \ \ 
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Figure 1. The final eccentricity as a function of strength of the 
wind asymmetry and the mass of the primary and secondary. The 
upper panel gives the case for equal-mass binaries, and the lower 
panel shows the results for secondaries that are much smaller 
than the primary (the results for even smaller secondaries look 
similar). The solid curves follow the component proportional to 
C± and the dashed curves follow the component in the direction 
of Lq X b. The final value of the eccentricity is given by e = 
|C|| -HAc^ oCi,o-hALoxb(iob)"^Lo xb| 
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Figure 2. The required velocity kick imparted to a single star 
required to disrupt an initially circular binary as a function of the 
masses of the two stars in the system. The curves are labeled with 
the mass of the primary. The required velocity kick increases as 
the mass of the secondary decreases because the contribution of 
the wind to the equation of motion of the binary is proportional 
to M2 from Eq. |[TJ. This is a result of the separation of the system 
into the motion of the centre of mass and the motion about the 
centre of mass. 



the resonance structure of multiple planet systems (see also 
iDebes fc SigurdssonllioO^ ) . possibly resulting in the ejection 
of smaller planets from the system. Finally the white dwarf 
kicks could play an important role in the binary evolution 
and especially double degenerate binaries, yielding systems 
that would otherwise be difficult to form. 



ACKNOWLEDGMENTS 

I would like to thank Harvey Richer and Saul Davis for use- 
ful discussions and Camp Byng where some of the work was 
completed. The Natural Sciences and Engineering Research 
Council of Canada, Canadian Foundation for Innovation and 
the British Columbia Knowledge Development Fund sup- 
ported this work. Correspondence and requests for materi- 
als should be addressed to heyl@phas.ubc.ca. This research 
has made use of NASA's Astrophysics Data System Biblio- 
graphic Services 



REFERENCES 

Davis S., Richer H. B., Coffey J., Anderson J., Brewer J., 
Fahlman G. G., Hansen B. M., Hurley J., Kalirai J. S., 
King I. R., Reitzel D., Rich R. M., Rich M. R., Shara 
M. M., 2006, in American Astronomical Society Meet- 
ing Abstracts Vol. 209 of American Astronomical Soci- 



